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We obtain a family of uniformly propagating hole-solutions to the complex Ginzburg-Landau 
equation which differ from the well-known Nozaki-Bekki holes. They describe the spatial orga- 
nization and velocity of the dominant local structures in the spatiotemporal intermittent regime. 
We discuss the phenomenology of these intermittent states in terms of the properties of the local 
structures. 
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A proper understanding of spatiotemporal chaos, i.e., 
deterministic chaos occurring in extended systems that 
are driven from equilibrium, is lacking. Since the num- 
ber of effective degrees of freedom diverges with the sys- 
tem size, most of the tools developed for low-dimensional 
systems are inapplicable. Moreover, these tools do not 
provide a proper framework to describe the spatial orga- 
nization of extended chaos. In many cases, the dynamical 
states appear to be built up from local, almost particle- 
like objects with well-defined dynamics and interactions 
(Fig. Ib-d). A description of spatiotemporal chaos in 
terms of these structures is therefore desirable 

In this Letter we will investigate the local structures 
that appear in the spatiotemporal intermittent regime of 
the complex Ginzburg-Landau equation (CGLE): 
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This amplitude equation describes pattern formation 
near a Hopf bifurcation and has been applied to de- 
scribe patterns occurring in, e.g., fluid-convection, Fara- 
day waves, optical systems, chemical oscillations, and 
turbulent flow past a wake |^ . As a function of the coeffi- 
cients ci and C3 , which are determined by the underlying 
physical problem, behavior ranging from completely reg- 
ular to strongly chaotic has been found (Fig. la). 

In the spatiotemporal intermittent regime |^,^ , a plain 
wave attractor coexists with a chaotic attractor; most 
initial conditions evolve to the latter (Figs. Ib-d). The 
typical states consist of patches of plane waves, separated 
by various "holes", i.e., local structures characterized by 
a depression of |^|. Similar intermittent states have been 
reported for the damped Kuramoto-Sivashinsky equa- 
tion, coupled cole's, Rayleigh-Benard convection and 
the printer instability [||^. It has been suggested that 
spatiotemporal intermittency should occur generally in 
the transition route from laminar to chaotic states, and 
the phenomenology suggests a relation to directed perco- 
lation 

The local structures in the intermittent regime can be 
divided into two groups depending on the wavenumbers 
qi and Qr of the asymptotic waves they connect. The 
quasi-stationary structures in Fig. Ic and Id have qi ^ qr 



and are related to the intensively studied Nozaki-Bekki 
holes However, the dominant local structures have 
velocities and asymptotic wavenumbers that are incom- 
patible with the Nozaki-Bekki holes jj]. For example, the 
holes shown in Fig. lb all have qi wg,. « since the plain 
waves (gray areas) have wavenumbers close to zero; sim- 
ilar holes are the dominant structures in the bi-chaotic 
regime. In the following we shall characterize these holes 
and their dynamical properties, and discuss their rele- 
vance for the chaotic states of the CGLE. 

In Fig. 2, \A\, the complex phase (arg(yl)) and the local 
wavenumber q:—dx arg(yl) for the left lower part of Fig. 
lb are shown. The wavenumbers of the laminar patches 
are quite close to zero, while the cores of the local struc- 
tures are characterized by a sharp phase-gradient (peak 

in q) and a dip of The holes propagate with a speed 
of 0.95 ± 0.1 and either their phase-gradient spreads out 
and the hole decays, or the phase-gradient steepens and 
the hole evolves to a phase-slip. As a first step in describ- 
ing these local structures which have a slowly evolving 
velocity and spatial structure, we will study the coher- 
ent structures, i.e., structures with fixed spatial structure 
and velocity. 

For the ID CGLE, coherent structures have been de- 
scribed in a simple framework fl^ . By substituting an 
ansatz for a uniformly propagating solution of the form 
A{x,t) = e-^'^^A{Cl into the CGLE ^ {i-^x~ vt), we 
obtain a set of coupled first order ordinary differential 
equations (ODE's) 
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[-1 - iuo + {I - ic.i)a^ - vz] , (2b) 



where a := \A\ and where the complex quantity z is de- 
fined as \n{A) —: n ~\- iq. Equation (^) is equivalent 
to two real valued equations, so (0) can be seen as a 3D 
real-valued dynamical system [ [13[ . Plain waves corre- 
spond to fixed points of (||), and the hole solutions we 
are interested in correspond to orbits connecting these 
fixed points. A rather complete study of the heteroclinic 
orbits which describe for instance the Nozaki-Bekki holes 
{qi ^ qr) has been made [Ol. We are here interested in 
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local structures that have qi= Qr, i.e., homoclinic orbits 
of the ODE'S (|). 

In general the ODE's have uj and v as free param- 
eters, but since the wavenumber in the laminar patches 
is approximately zero, we demand qi—qr — 0, which fixes 
UJ ~ ~C3. The fixed point at (a, z) = (1,0) corresponds 
to the q — plain waves and has a ID outgoing manifold 
and a 2D spiraling ingoing manifold (Fig. 3a). To create 
a homoclinic orbit, we have to connect these manifolds. 
This amounts to satisfying only a single condition, and 
since we have one free parameter (w), we can expect a 
discrete set of homoclons for qi — g,. — 0. Performing a 
simple numerical integration of (H) and adjusting the free 
parameter v, we obtain a homoclinic orbit for v « 0.916; 
the corresponding coherent structure (Fig. 3b) will be 
referred to as a homoclon. For qi = qr ^ 0, i.e., uj ^ — C3 
one can obtain similar homoclinic orbits, so in fact there 
exists a one-parameter family of homoclons. 

In Fig. 3b the homoclons are compared to the local 
structures in the intermittent regime. The slight devia- 
tion between this particular local structure and the ho- 
moclon is mainly due to the fact that the local structure 
has a slowly evolving shape. The longer the lifetime of the 
holes, the better the fit is to the homoclons. For nearby 
values of ci and C3 one finds similar correspondences be- 
tween the q = homoclons and the local structures. 

It is instructive to compare the homoclons with the 
Nozaki-Bckki holes [|l2|. The Nozaki-Bekki holes con- 
tradict naive counting arguments and have been shown 
to be structurally unstable, while the homoclons satisfy 
the counting arguments, and are structurally stable. One 
can verify that for ci = 0.6, 03 = 1.4, a Nozaki-Bekki hole 
with qi~Q has qr « 0.837 and a velocity of 1.673, which 
is completely different from the local structures and the 
homoclons. In some regimes, in particular for negative 
ci, the Nozaki-Bekki holes are dynamically relevant, and 
as they are sources for waves with q ^ 0, we obtain grain 
boundaries between q =^ and g « waves (Fig. lc,d). 
In the limit where qi = qr, the width of the Nozaki-Bekki 
holes diverges, so we can only conclude that there are 
two distinct types of hole solutions: heteroclinic Nozaki- 
Bekki holes and homoclons. 

Homoclinic orbits of (^) describe at least two other 
cases. Since the homoclinic orbit occurs as a saddle-focus 
loop, there are multi-circuit loops and a limit-cycle for 
nearby values of t; or w fl^ . The limit cycle collapses via 
a Hopf bifurcation on another fixed point of the ODE's 
(^ for sufhciently small v. The (quasi)periodic states 
of the CGLE that have been called "compression waves" 
I p^ are described by these limit-cycles. Furthermore, 
the "riding chaos" recently reported in has been de- 
scribed in the context of homoclinic orbits for the phase 
equation of the CGLE, and these orbits also occur in 
the ODE's (^. For example, the uniformly propagating 
structures shown in Fig. II of occur in a background 
of wavenumber 0.038 for cl = 1.75, c3 = 0.8. 



The ODE's for the full CGLE (|) admit a homoclinic 
orbit for vRiO.557, that describes the riding chaotic states 
perfectly. For the same asymptotic wavenumbers and 
values of the coefficients, there exists also a homoclinic 
orbit for i;«I.I24; this orbit is similar to the homoclons 
and cannot be obtained from the phase equations. More 
detailed studies of the homoclinic orbits are underway 

EMI-. 

We will now return to the homoclons which appear in 
the spatiotemporal intermittent regime and study their 
stability. Following an isolated, slightly perturbed homo- 
clon, one immediately finds that they are weakly unstable 
. In a similar fashion to the local structures found in 
the intermittent state, they either slowly decay or grow 
out to a phase-slip (Fig 4a-d) . In the phase-space of the 
CGLE (0) we can think of the homoclons as unstable 
equilibria separating plain waves and phase-slips. 

In Fig. 4a we show the decay of a slightly perturbed 
homoclon (perturbation of order 10~^). The total phase- 
difference A6 across the decaying homoclon is conserved 
since there are no phase-slips; for ci — 0.6 and C3 = 1.4, 
A9 w 3.24. During the decay, the wavenumber-peak, 
amplitude-dip and apparent velocity decrease; for long 
times, the dynamics crosses over to a slow phase-diffusion 
by which A9 is smeared out. In Fig. 4b, the evolution 
of a slightly perturbed homoclon towards several phase- 
slips is shown. The wavenumber-peak and amplitude- 
dip slowly grow, and at f « 117 the first phase-slip oc- 
curs, from which a typical spatiotemporal intermittent 
state nucleates. Before this phase-slip, the wavenumber 
acquires a negative peak in order to conserve the total 
phase-difference across the structure (Fig. 4c). Both 
these peaks diverge at the phase-slip event, and just af- 
ter the phase-slip, the winding-number J dxq/2n has de- 
creased by 1. Therefore, the negative phase-bump that 
corresponds to the new left moving hole, is quite steep 
(Fig. 4d), and this hole will quickly grow out to a new 
phase-slip, from which a strong right moving hole is gen- 
erated etc. When we quench ci and C3 in the direction 
of the transition to plain waves, these zigzag motions of 
the holes become very dominant (Fig. 5a). 

The group- velocity for the plain waves is close to zero. 
Therefore, in the co-moving frame of a right moving hole, 
the wave to the right has its group-velocity pointing in- 
ward and the wave to the left has it pointing outward; 
hence the homoclons are neither sources nor sinks. Sup- 
pose we have a homoclon with a positive phase-difference 
A9, moving to the right into a plain wave with wavenum- 
ber qr- Because they are unstable, their dynamical fate 
depends strongly on qr- When qr > 0, this leads to 
the "winding up" of the homoclon (increase of A6), and 
a phase-slip is generated in which the total winding- 
number is decreased. The larger qr is, the faster this 
process goes, and as a result, the total number of phase- 
windings in the system is driven to approximately zero. 
When qr < 0, or equivalently, when a left moving ho- 
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moclon invades a state with positive g, this leads to 
the "winding down" of the homoclon, which then de- 
cays. Its associated is smeared out and this drives 
the wavenumbers in the laminar patches to zero. The 
winding up and down is shown in Fig. 5b, where we 
follow the evolution of an isolated homoclon in a back- 
ground state with wavenumber 0.05. This strong sen- 
sitivity to the asymptotic waves has as a consequence 
that the fluctuations in the laminar patches of the inter- 
mittent regime have a strong effect on the homoclons. 
Small wave-packets, resulting from decaying homoclons, 
can trigger the decay or phase-slip of the propagating 
holes; this might result in correlations in the intermit- 
tent state over long times and distances. 

The homoclons have only one weakly unstable eigen- 
mode. Consider the example illustrated in Fig. 6a. 
The initial condition here has |^| = 1 and a triangular 
wavenumber profile. In general, there is no reason to ex- 
pect such an initial condition to evolve to a homoclon, 
but by adjusting only one parameter in the initial con- 
dition, this is precisely what happens; for the example 
shown in Fig. 6a, the height of the triangle was set to 
0.437754 while its width was 20. Increasing the height 
or width of the initial wavenumber-blob leads to a steep- 
ening hole like in Fig. 4b, while a decrease leads to a 
decaying hole like in Fig. 4a. The total phase-difference 

across the triangular initial condition is 8.76, so a big 
positive wavenumber-packet is emitted from the homo- 
clon. We have also tried this procedure for other initial 
conditions, and found that in general the adjustment of 
only one free parameter in the initial conditions leads to 
the generation of a long living homoclon; the dimension 
of the unstable manifold is therefore one. The contrac- 
tion towards this unstable manifold is rather fast. This 
is reflected in the dynamics in the intermittent regime; 
the local structures evolve to a reasonable approximation 
along the ID unstable manifold of the homoclons. 

To substantiate this claim, we considered the relation 
between the values of the extrema of q and the corre- 
sponding local minima of |^|. As shown in Fig. 6b, these 
quantities are strongly correlated. This indicates that a 
one-parameter family of profiles of A is dominant. Note 
the increase of the density around values of a and q which 
correspond to the coherent homoclon; being an equilib- 
rium state, the dynamics spend a relative long time there. 

In Fig. 6c, d we collapsed the profiles which correspond 
to certain isolated minima of |^|. The profiles shown in 
Fig. 6c correspond to the coherent homoclon, while the 
profiles in Fig. 6d correspond to states evolving towards 
a phase-slip. 

This suggests a phenomenological model in terms of 
moving "homoclon" particles possessing an internal de- 
gree of freedom that parameterizes the location on the 
unstable manifold. More detailed studies of this are un- 
derway pof . 



In conclusion, we have described a new class of coher- 
ent solutions which occur in the spatiotemporal intermit- 
tent regime and which are intimately connected to phase- 
slips; they are therefore one of the prime local structures 
of the ID CGLE. 
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FIG. 1. (a) "Phase-diagram" of the CGLE. For small ci and C3 all initial conditions evolve to plain waves. In the 
intermittent regime, a plain wave attractor and a chaotic attractor coexist. Beyond the full curve ciC3 = l, all plain waves are 

linearly unstable and all states are spatiotemporal chaotic. At a zero of A the complex phase is undefined and phase-slips occur 
(see Fig. 2); the chaotic state is then called defect-chaos. When A has no zeroes we speak of phase-chaos. In the bi-chaotic 
regime, a defect and phase-chaotic attractor coexist [4,11]. (b)-(d) Space-time plots (over a range of 200 x 150) of |^| (black 
corresponds to |^| = 0) showing chaotic states in the spatiotemporal intermittent regime, for coefficients (01,03) = (0.6, 1.4) (b), 
(0,1.8) (c) and (0,1.4) (d). 



FIG. 2. Space-time (60x50) plots of the left-lower part of Fig. lb, showing \A\, arg(A) and q:=dx arg(A) in detail. 

FIG. 3. (a): The homoclinic orbit of the ODE's (2) in a,q,K, space, (b) The amplitude and wavenumber profile of the 

corresponding coherent structures (curves) . The circles correspond to a local structure obtained from simulations of the CGLE 
in the spatiotemporal intermittent regime. The bump in a to the right of the core corresponds to the spiraling motion on the 
incoming manifold of the fixed point (o, q, k) = (1, 0, 0). 

FIG. 4. (a,b) Evolution of the wavenumber profiles of slightly perturbed homoclons. Consecutive time-slices have a time 
difference of 5. (c,d) Wavenumber profiles (full curves) and winding-number J dxq/2-K (dashed curves) just before (c) and after 
(d) the first phase-slip. 

FIG. 5. (a) The dominance of zigzagging holes near the transition to plain waves ci = 0.6, C3 = 1.2 for space x time 
= 512 X 1000. (b) The evolution of a homoclon in a background state with wavenumber 0.05, for ci = 0.6, C3 = 1.4 and space x 
time = 512 x 250. 



FIG. 6. (a) The waveuuuibcr profiles of a "triangular" initial condition evolving to a homoclon. (b) Collapse of the extrema 
of a and q in the intermittent regime, (c) Amplitude profiles of j^l with minima of |^| around the homoclon value, (d) Similar 
profiles for steeper minima. 
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